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ABSTRACT 


A  mathematical  computer  model,  "HARRIER",  of  the  AV-8A  Harrier  air¬ 
craft  has  been  developed  utilizing  the  concept  of  quasi-stationary  flight 
in  order  to  determine  quasi-equilibrium  values  of  flight  variables  and 
controls  during  accelerating-decelerating  flight.  Using  available  aero¬ 
dynamic  data  via  look-up  schemes,  "trim-maps"  are  computed  which  show  the 
values  of  the  longitudinal  flight  and  control  variables  as  functions  of 
airspeed  and  thrust  angle.  These  maps  have  been  shown  to  be  useful  tools 
in  piecing  together  transition  profiles,  taking  in  consideration  various 
requirements  and  constraints. 

The  final  decelerating-descending  phase  of  the  present  U SMC  approach 
profile  has  been  analyzed  in  detail  with  the  aid  of  the  trim-maps.  The 
technique  of  flying  the  transition  as  suggested  by  the  trim-maps  has  been 
found  in  good  agreement  with  the  techniques  actually  used  under  the  given 
constraints.  This  established  confidence  in  both  the  mathematical  model 
and  in  the  trim-maps  as  tools  for  studying  alternate  transition  profiles. 
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CHAPTER  I 
INTRODUCTION 

Over  the  past  decade  Vertical/Short  Takeoff  and  Landing  (V/STOL)  aircraft 
have  been  of  increasing  interest  to  the  aerospace  community.  One  of  the  most 
critical  phases  of  V/STOL  flight  is  the  transition  from  wing  borne  to  jet  borne 
flight.  The  stability  of  a  V/STOL  aircraft  normally  changes  significantly  in 
transition  and  becomes  unstable  at  hover.  This  creates  a  high  pilot  workload 
and  therefore  a  potentially  dangerous  situation.  Therefore,  investigating 
transition  techniques  and  flight  paths  is  of  great  importance  when  studying 
V/STOL  aircraft. 

The  purpose  of  this  report  is  to  develop  and  demonstrate  a  mathematical 
method  for  studying  the  transition  performance  and  control  of  highly  non-linear 
V/STOL  aircraft.  Although  the  sample  aircraft  chosen  is  the  AV-8A  Harrier, 
the  analysis  used  is  valid  for  any  V/STOL  aircraft.  The  purpose  of  this  work 
is  to  demonstrate  the  usefulness  of  the  method  to  the  analysis  of  the  Harrier 
and  to  provide  some  general  guidelines  on  Harrier  transition  techniques. 

The  method  is  based  on  a  non-linear  longitudinal  force  and  moment  balance 
model  of  the  Harrier  aircraft  in  transition  flight.  The  concept  of  a  force 
balance  model  is  extended  to  include  the  inertial  forces  produced  by  constant 
acceleration  along  a  prescribed  flight  path.  This  model  represents  "snapshots" 
of  the  force  and  moment  balance  while  accelerating  along  the  flight  path.  The 
concept  of  "trim"  is  expanded  to  mean  the  combination  of  controls,  at  a  given 
flight  condition  (altitude,  weight,  c.g.  position  and  velocity),  which  is  needed 
to  establish  force  and  moment  balance  including  the  inertial  components  due  to 
acceleration.  If  it  is  desired  to  examine  a  trim  point  where  no  flight  path 
acceleration  is  taking  place  then  "trim"  assumes  its  conventional  meaning  of 
an  aircraft  in  steady  unaccelerated  flight.  This  model  is  used  to  create  a  "trim- 
map"  of  the  variables  of  interest  in  the  transition  speed  range.  With  such 
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"trim-maps"  the  performance  and  control  limitations  of  the  aircraft  may  be 
studied  in  a  logical,  organized  fashion.  The  trim-maps  can  be  used  to  identify 
candidates  for  the  best  transition  flight  path  and  conversion  program,  taking 
into  consideration  many  factors,  such  as  fuel  consumption,  time  duration,  pilot 
workload  and  visibility,  control  margins,  and  engine  operating  limits. 

Mapping  the  trim  points  is  a  convenient  way  to  present  the  result  of  the 
force  balance  model  because  in  V/STOL  aircraft  there  is  generally  a  redundancy 
in  controls  during  conversion/transition. 

In  the  Harrier,  the  pilot  has  four  longitudinal  controls  available: 

1.  Longitudinal  stick 

2.  Throttle 

3.  Flaps 

4.  Jet  deflection  angle 

Figure  1  illustrates  how  different  combinations  of  longitudinal  stick, 
(controlling  the  attitude),  jet  deflection  angle  (nozzle  position),  and 
throttle  position  may  combine  in  different  ways  to  create  the  same  flight 
conditions.  With  the  flaps  down,  the  Harrier  can  fly  straight  and  level  at 
220  ft/s  at  any  angle  of  attack  (c)  between  15.5°  and  -5.0°  with  the  nozzle 
deflection  (9^)  varying  from  0°  to  80°.  Varying  the  flap  setting  (6^)  would 
further  increase  the  variations  in  a  and  0^ .  The  goal  of  this  study  is  to 
establish  the  relationships  among  the  longitudinal  controls,  including  control 
limitations,  so  that  transition  flight  paths  and  piloting  techniques  may  be 
compared  as  to  their  relative  merits. 

First,  in  the  rest  of  this  Chapter,  some  background  and  a  brief  description 
of  the  AV-8A  Harrier  are  given.  Chapter  2  discusses  the  development  of  the 
mathematical  model.  Finally,  in  Chapter  3  the  resulting  trim-maps  are  discussed 
and  an  example  transition  flight  path  is  given  to  illustrate  how  the  trim-maps 
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are  utilized.  Henceforth,  HARRIER  in  capital  letters  refers  to  the  computer 
program  . 

Background 

In  studying  the  aforementioned  problem,  much  effort  was  expended  in  re¬ 
searching  methods  to  theoretically  determine  V/STOL  transition  flight  paths 
(References  12  through  15).  Although  many  papers  have  been  written  on  this 
general  subject,  it  was  found  that  none  of  the  studies  accounted  for  an  air¬ 
frame  that  had  severe  control  power  limitations  and  coupling  between  engine 
thrust  and  aerodynamic  lift.  The  fact  that  the  lift  force  is  a  function  of 
gross  thrust  and  thrust  deflection  angle  in  addition  to  the  angle  of  attack 
vastly  increases  the  analytical  problem.  In  addition,  when  the  aircraft  has 
marginal  control  power  within  a  certain  speed  range,  the  problem  of  finding 
the  best  way  to  transition  the  aircraft  from  forward  to  hovering  flight  is 
compounded.  This  is  because  the  performance  and  the  stability  and  control 
characteristics  must  be  considered  simultaneously.  The  path  for  a  minimum 
power,  minimum  fuel  consumption,  or  minimum  time  approach  may  be  (and  is  for 
the  AV-8A)  severely  constrained  by  the  longitudinal  control  power  available. 
Closed  form  flight  path  optimization  becomes  very  complicated  when  the 
above  mentioned  problems  are  included. 

Another  problem  with  classical  optimization  in  general  is  that  once  the 
optimal  flight  path  is  found,  only  the  pilot  technique  required  to  fly  this 
path  is  studied.  It  becomes  difficult  to  research  the  "optimum  pilot  tech¬ 
nique"  in  any  sense  in  this  manner. 

The  trim-maps  derived  and  discussed  in  this  report  incorporate  all  the 
significant  features  of  the  aircraft.  They  are  equally  applicable  to  gen¬ 
eral  theoretical  aircraft  or  to  any  particular  V/STOL  aircraft.  The  only 


drawback  of  the  trim-maps  is  that  they  do  not  directly  show  how  to  get  from 
one  point  to  another.  In  most  cases  this  is  relatively  straightforward. 
They  do,  however,  vividly  show  minimum  fuel  consumption  (thrust)  approaches 
and  pilot  technique  required.  Control  margins  and  effectiveness  are  also 
indicated . 

The  interpretation  of  the  calculated  trim-maps  for  the  AV-8A  is  con¬ 
tained  in  Chapter  3. 


Description  of  the  AV-8A  Harrier 

The  AV-8A  Harrier  is  a  subsonic,  single  cockpit,  single  engine  turbofan 
aircraft  with  V/STOL  capabilities  which  arise  from  being  able  to  rotate  the 
thrust  vector  from  0°  (thrust  aligned  with  the  x  axis)  to  98.5°  by  means 

O 

of  four  nozzles,  two  on  each  side  of  the  aircraft.  The  two  front  nozzles 
exhaust  the  low  pressure  bypass  air  and  the  two  rear  nozzles  exhaust  the  hot 
turbine  air.  Gross  thrust  is  distributed  almost  equally  between  the  forward 
and  aft  nozzles  and  the  resultant  thrust  vector  acts  close  to  the  eg. 

The  Harrier  has  a  tapered  swept  wing  mounted  on  top  of  the  fuselage 
with  negative  geometric  dihedral.  Similarly,  the  horizontal  tail  too  is  ta¬ 
pered,  swept  and  has  negative  geometric  dihedral.  There  is  a  single  ventral 
fin  mounted  below  the  vertical  stabilizer  at  the  tail.  The  aircraft  has 
half-span  trailing  edge  flaps  which  are  normally  deflected  50°  for  takeoff 
and  landing  and  0°  (flaps  up)  for  cruise.  There  is  also  a  speedbrake,  mount¬ 
ed  hinge  forward,  beneath  the  tail  and  just  forward  of  the  ventral  fin.  Max¬ 
imum  speedbrake  deflection  is  66°.  Nominal  speedbrake  positions  are  0°  (speed- 
brake  stowed)  for  cruise  and  2S°  for  takeoff  and  landing  configuration  (un¬ 
dercarriage  down) .  With  the  undercarriage  down  the  speedbrake  cannot  be 
moved  from  the  25°  position. 


The  AV-8A  is  approximately  45  ft  long  and  has  a  25-ft  wingspan.  Empty 
operating  weight  is  approximately  12000  lbs.  Maximum  takeoff  gross  weight 
is  24600  lbs  and  design  landing  weight  is  16900  lbs.  Maximum  vertical  take¬ 
off  gross  weight  is  175 00  lbs.  Maximum  gross  engine  thrust  is  slightly  less 
than  20000  lbs. 

Figure  2  shows  a  three  view  drawing  of  the  AV-8A  Harrier. 

One  aircraft  configuration  was  examined  in  detail  in  this  study; 

a)  flaps  down  50° 

b)  undercarriage  down 

c)  speedbrake  extended  30° 

As  previously  mentioned,  actual  speedbrake  deflection  in  takeoff  and 
landing  is  25°.  However,  because  the  wind  tunnel  data  were  available  only 
for  30°  deflection,  the  model  assumes  a  30°  speedbrake  deflection. 

Most  of  the  data  used  in  modelling  the  AV-8A  is  contained  in  five 
McDonnell  Douglas  Corporation  reports.  References  16  through  20.  Reference 
16  contains  the  vast  majority  of  the  aerodynamic  data  on  the  Harrier. 
Reference  17  contains  reaction  control  system  data,  while  Reference  18 
details  the  jet-induced  lift  loss  effects.  The  model  of  the  reaction  con¬ 
trol  system  described  in  this  report  was  derived  based  on  Nave's 
report  (Ref.  7.)  which  contains  the  description  of  a  simulation  of  the 
dynamics  of  the  AV-8A  in  transition. 
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CHAPTER  II 


AV-8A  HARRIER  MODEL  DEVELOPMENT 
Quasi-Equilibrium  Flight  Model 

For  the  purposes  of  producing  a  set  of  trimmaps  a  quasi-equilibrium  flight 
model  is  needed.  In  this  study  equilibrium  flight  is  defined  as  steady  un¬ 
accelerated  flight  along  a  straight  path  (not  necessarily  level).  This 
means: 

u  ,  v  ,  and  w  are  constant 

V  V  V 

uv=vv=wv=p=q=r=p=q=r=0 
Quasi- equilibrium  flight  is  defined  similar  to  equilibrium  flight  but  allows 
for  constant  linear  accelerations.  Having  both  velocities  and  accelera¬ 
tions  constant  in  the  force  and  moment  balance  equations  means  that  only  a 
particular  instant  of  time,  or  a  "snapshot”  of  am  aircraft  accelerating 
along  a  flight  path,  is  examined.  The  instant  of  the  snapshot  is  the  time 
at  which  the  aircraft  accelerates  through  the  specified  velocity.  In  the 
quasi-equilibrium  case; 

u  ,  v  ,  w  ,  u  ,  v  ,  w  are  constant 
v  v  v  v  v  V 

P  »  q  >  r  ,  p  ,  q  ,  r  are  small  and  assumed  zero 
These  assumptions  simplify  modelling  the  aircraft  by  eliminating  time  de¬ 
pendence,  so  only  the  algebraic  equations  of  force  and  moment  equilibria 
need  be  solved. 

General  equations  derived  with  the  assumptions  above  are  presented  in 
Appendix  A.  Because  this  study  only  considers  flight  in  a  vertical  plane, 

I  *  $  *  5  3  yj*  0; 
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rife 


r 

f 


furthermore  the  ambient  wind  is  considered  to  be  zero: 

6  '  vB  ■  0 

Using  these  assumptions,  the  seven  equations  required  to  model  the  Harrier 
in  transition  are  Equations  A-3,  A-5,  A-13,  and  A-16  from  Appendix  A. 


u„  =  (1/W) (-Dcosa  +  Lsina  +  Tcos0.)  -  sin0 
B  3 

w  ■  (1/W)  (-Dsina  -  Lcosa  -  Tsin0 .  -  TQ  -  Ta  )  +  cos0 
B  2  6A  9F 

L  =  (1/2)PV*SCL 

D  =  Cl/2)PV2SCn  ♦  D 

*T  *  u  cos0  ♦  w  sin8 
IB  B 

i-  »  -u  sin8  +  w  cos0 

IB  d 

M  =  ((l/2)pvJscC_  +  Tq  x0  +  Tfi  xfl  +  (T.  x.  +  T.  x.  )sin0. 

w  ,  0  4  jp  Jp  J  ^  J  ^  J 


CD 

(2) 

(3) 

(4) 

(5) 

(6) 


F  F 


A  A 


+  x  Lcosa  +  x 
mr  -mx 

Dsina) /I 

yy 

(7) 

i  M  =  0  in  quasi-equilibrium  flight 

a  =  0  -  yc  (with  no 

ambient 

winds ) 

(8) 

T  =  (V./VJ2qSj 

(9) 

C  *  f(a,  ~-l—  , 

1  <s 

Vj_ 

v« 

9i  -5f  -V 

(10) 

CD  '  -7-  ■ 

pv;s 

V. 

' 

Voo 

9j  -  V 

(11) 

D  =  f(V  ,T) 

ram  00 

C  a  f (a,  — \ —  , 

" 

V. 

/  * 

00 

9j  >SF  -V 

(12) 

(13) 

T  ’  «T.  V.,  V 
A 

(14) 

T0R  *  f<T>  V 

F 

(15) 

T  -  f  (T,  VJ 

JA 

(16) 

T  a  f  (T,  VJ 

(17) 
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and  xa  ,  xn  ,  x.  ,  x  ,  x  are  all  moment  arm  lengths  and  functions 
0F  ®A  JF  JA  mr 

of  the  longitudinal  eg  position.  They  are  assumed  to  be  constant  for  each 
flight  profile.  S,  W,  I  ,  and  c  are  wing  area,  aircraft  gross  weight, 
pitching  moment  of  inertia,  and  mean  aerodynamic  chord,  respectively;  p  is 
the  atmospheric  density  and  is  held  constant  at  its  sea  level  value. 

Equations  (1)  and  (2)  define  the  body-axis  accelerations  of  the  air¬ 
craft.  Equations  (3)  and  (4)  compute  the  aerodynamic  lift  and  drag  of  the 
aircraft.  Equations  (5)  and  (6)  comprise  the  axis  transformation  from  the 
body-axis  system  to  the  local  inertial  axis  system.  Finally,  Equation  (7) 
is  the  equation  for  pitch  acceleration.  The  actual  functions  C^,  C^,  and 
are  quite  complex  and  shall  be  discussed  later  in  this  chapter. 


Creation  of  the  Trim-Maps 


Equations  (1)  through  (7),  with  the  relations  expressed  by  Equations 

(8)  through  (17),  allow  for  the  solution  of  7  unknowns.  However,  a  total 

of  25  unknown  parameters  and  variables  appear  in  the  17  equations  after  the 

aircraft  physical  characteristics  (S,  c,  W,  I  ,  eg,  V*  X0  *  V  Xjr  ' 

A  F  A  F 

are  given.  These  unknowns  are: 

V  V  V  V  L’  D’  a’  9’  T’  T0a’  T6f’  CL’  V  Dram'  Cm'  TjF’ 

TjA*  v  v  v  v°°’  yv**  y*  v  li 

Therefore,  after  the  physical  characteristics  are  set,  8  more  parameters 
or  variables  must  be  specified  or  determined. 


For  each  trim-map  the  flap  position,  the  vertical  flight  path  angle, 

Y  ,  and  the  accelerations  of  the  aircraft  in  the  flight  path  (velocity)  axes, 

u  and  w  (the  subscript  c  indicates  commanded  quantities)  are  specified 
c  c 


and  held  constant.  Airspeed,  V  ,  and  jet  nozzle  deflection  angle,  0.,  are 


varied  to  map  the  speed  range  where  quasi-equilibrium  flight  may  occur  at 
each  0j.  These  six  variables  and  parameters  are  inputs  to  the  HARRIER 
computations . 

The  body  axis  accelerations  in  Equations  (1)  and  (2)  are  substituted 
by  flight  path  axis  acceleration  inputs  via  the  following  axis  transformation 


•  u  cosa  -  w  sma 

V  V 

c  c 

(from  A- 13) 

(18) 

=  u  sina  +  w  cosa 

v  V 

(19) 

c  c 


This  substitution  does  not  affect  the  number  of  unknowns. 

Two  more  variables  are  actually  inputs  as  determined  by  the  following 
two  equations: 

*1  =  cosYc  (20) 

(from  A-8) 

fi  =  -v~  sinYc  (2D 

C 

Substituting  Equations  (18),  (19),  (20),  and  (21)  into  the  left  sides 

of  Equations  (1),  (2),  (5),  and  (6)  respectively  and  rearranging  Equations 

(1)  through  (7) ,  with  q  =  (l/2)pV^  , 

0  =  (1/W)  (-Dcosa  +  Lsina  -  Tcos  0. )  -  sin  9  -  (u  /g)cosa  +  (w  /g)sina 

^  vc  vc 

(22) 


0  =  (1/W)(-Dsina  -  Lcosa  -  Tsin8.  -  T.  -  TQ  )  +  cos0 

3’  \  0F 

-  (uy  / g)sina  -  (wy  /g)cosa 
c  vc 


0  *  L  -  qSCL 

0  *  D  -  qSCn  -  D 

n  D  ram 


(23) 

(24) 

(25) 
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(26) 


0  =  -uDcos  9  -  wnsin9  +  V  cosy 
B  B  oo  !c 

c 


0  =  u  sin  9  -  wncos9  -  V  siny 

D  D  00  C 

c 

0  =  (qScC  +  T_  xQ  Tq  xa  +  (T.  x.  +  T.  x.  )  sin9 . 

"  9f  6f  9 a  9a  h  h  >A  >A  9 


(27) 


+  x  Lcosa  +  x  Dsina)/I 
mr  mr  yy 


(28) 


For  determining 


a,  C.  ,  CL,  D  ,  C  ,  V./V  ,  TQ  ,  T„  ,  T.  ,  and  T. 

L’  D*  ram’  m*  ]'  “»*  9^  8p’  jp’  j A 

the  relationships  expressed  in  Equations  (8)  through  (17)  are  used,  utiliz¬ 
ing  look-up  schemes  and  interpolation  where  necessary,  making  a  total  of  17 
independent  relationships  available  to  calculate  17  variables.  The  quasi¬ 
equilibrium  "trim"  points  for  each  variable  of  interest  can  be  calculated 
and  plotted  as  functions  of  Vm  and  0^.  These  plots,  one  for  each  variable 
(L/W,  D/W,  T/W,  9,  6U,  u  ,  w  )  are  called  trim-maps. 

Appendix  B  describes  the  method  used  to  solve  the  system  of  equations. 


Aerodynamics 

Modelling  the  aerodynamics  of  the  Harrier  aircraft  is  a  complex  task 
because  of  wing/jet  aerodynamic  interaction.  This  interaction  is  a  result 
of  the  air  flow  at  the  engine  inlet  and  of  the  jet  exhaust  from  the  main 
nozzles  interfering  with  the  freestream  air  flow  at  the  wing.  The  inter¬ 
action  varies  strongly  with  the  jet  velocity  ratio,  V^/V^,  and  0_. ,  but  only 
slightly  with  a.  Essentially  all  aerodynamic  coefficients  except  speed- 
brake  effects  are  influenced  by  this  interaction.  The  effects  on  the  lift 

are  the  most  significant,  with  increasing  V./V  and/or  9.  decreasing  the 

J  00  J 

lift  produced  (for  constant  a).  This  results  in  a  significant  "lift-loss" 
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occurring  in  transition  from  forward  flight  in  the  airplane  configuration 

to  hover.  As  would  be  expected,  drag  generally  decreases  as  lift  decreases 

because  of  the  reduction  in  induced  drag.  As  V./V  and  0.  increase,  C 

s  J  00  1  m 

changes  from  negative  to  positive  and  C  becomes  neutrally  stable  or  unsta- 
ble.  In  the  majority  of  the  cases  these  terms  are  quite  non-linear  func¬ 
tions  of  V./V  and  9.. 

J  00  J 

The  interaction  effects  must  be  included  when  modelling  the  Harrier. 
Theoretical  models  representing  these  interactions  were  thoroughly  investi¬ 
gated  (References  1  through  4,  among  many)  only  to  find  that  because  of  the 
difficulty  and  poor  results  in  predicting  interaction  effects  all  sources 
recommend  the  use  of  experimental  data  when  available.  Therefore,  the  aero¬ 
dynamic  data  for  the  Harrier  which  include  interaction  effects  were  not 
modelled  analytically. 

The  basic  aerodynamic  data  are  presented  in  Reference  5  and  are  in  the 
form  of  carpet  plots.  Since  the  data  were  so  non-linear  and  functions  of 
three  variables  it  was  found  impractical  to  obtain  functional  relations  by 
means  of  curve  fitting.  Therefore,  the  aerodynamic  data  and  most  other 
functional  relationships  were  inserted  into  the  HARRIER  computer  program  in 
the  form  of  table  look-ups.  The  table  look-up  subroutines  give  the  value 
of  a  parameter  at  any  point  and  also  its  partial  derivative  with  respect 
to  each  of  the  inputs.  These  partials  are  needed  for  the  numerical  solution. 
Between  data  points  linear  interpolation  is  used  for  obtaining  the  value  of 
the  function  so  that  the  partial  derivatives  are  piecewise  constant  func¬ 
tions  of  the  arguments. 

The  components  of  C^»  and  used  to  model  the  Harrier  are: 
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c  *  c 

wing-body 


(a,  6.,  V./V.,  fiF)  ♦  C  (V./V,,  9.)A6h 

5h 


+  C.  (a,  5  )  +  c  (a,  V  ,  RPM) 

Lspeedbrake  inlet 


CD  *  CD  .  ,  ,  (0t’  ej»  Vj/V«»  V  +  CD  ...  (-a’ 

wing-body  J  J  speedbrake 

Ram  drag  is  computed  separately  and  is  therefore  not  included  in  C 


C  was  found  insignificant  and  is  therefore  neglected. 


H 


C  =  C 

m  m  .  ,  , 

wmg-body 


(a,  9.,  V./V*,  «p)  ♦  Cm  CV./V  a,  9.)A6h 

5h 


+  c  (a,  +  C 

mspeedbrake  “inlet 


(a,  V.,  RPM) 


where  =  6^  -  6.0° 

References  5,  6,  and  7  contain  the  data  which  determine  the  values  of  the 
above  coefficients. 


Wing-body  Aerodynamics 

The  wing-body  contribution  to  C  ,  Cn  and  C  with  flaps  down  is  present 

L  L)  IQ 

ed  in  Reference  5  in  the  form  of  carpet  plots  showing  Cr  as  a  function  of 

L* 

a,  and  C_  and  C  as  functions  of  C,  ,  for  various  V./V  and  9..  These  data 
’Dm  L  J  00  J 

were  taken  from  wind  tunnel  tests  of  a  powered  model  and,  therefore,  wing/ 
jet  interaction  effects  are  included.  Engine  thrust  and  moment  are  not  in¬ 
cluded  in  these  data  because  the  airflow  source  simulating  the  engine  was 
supported  separately  from  the  rest  of  the  model.  For  insertion  into  the 
HARRIER  computer  program,  CD  and  Cm  were  transformed  to  be  functions  of  a 

instead  of  C  .  Figures  3  through  14  display  these  data. 

L 

In  Reference  5  no  data  are  given  for  0  equal  to  zero  (airplane  con¬ 
figuration)  at  any  V./V  other  than  zero  (unpowered  model).  Therefore, 

J  00 
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nothing  is  known  about  the  wing/ jet  interaction  effects  for  0^=0.  To  in¬ 
vestigate  this  case,  the  carpet  plots  of  CL  ,  and 

wing-body  wing-body 

C  were  replotted  as  functions  of  9 .  for  a  constant  a  and  various 

wing-body  J 

Vj/V^.  Figures  IS  through  32  show  these  plots. 

From  these  plots  it  can  be  seen  that  for  constant  V  /V^  and  regardless 

of  angle  of  attack  and  approach  their  respective 

wing-body  u wing-body 

values  for  V./V^  *  0  as  9.  goes  to  zero.  This  indicates  that  for  C 

^  J  Lwing-body 

and  CD  there  are  negligible  wing/ jet  interactions  at  0 .  =0,  regard¬ 
wing-body  J 

less  of  V./V  . 

3  * 


While  CL  and  C, 

wing-body 

C  is  a  non-linear  function  of  9 . .  It  is  therefore  more  difficult 


-p  vary  almost  linearly  with  0 . , 

wing-body  ^ 


wing-body 
to  extrapolate  C 


to  9.  =  0.  However, 
mwing-body  ^  wing-body 


and 


"D  . 


clearly  show  that  wing/ jet  interactions  become  negligible  as  9. 


wing-body 

goes  to  zero  so  therefore  it  will  be  assumed  this  holds  true  for  C 
as  well. 


mwing-body 


As  a  result,  for  0^  =  0,  it  is  assumed  that  there  are  no  wing/ jet  inter¬ 
action  effects  regardless  of  the  value  of  V./V^.  Therefore,  in  Figures  3, 

7,  and  11  the  single  curve  represents  all  values  of  9^  at  V_. /V^  =0.  In  the 
rest  of  Figures  3-14  the  same  curves  as  on  Figures  3,  7,  and  11,  respective¬ 
ly, represent  the  values  of  ,  and  for 

wing-body  wing -body  wing-body 

0.  ■  0  for  all  V./V  . 

3  3  00 

Some  discrepancy  was  found  in  the  data  in  Reference  5  when  the  data 
from  the  aforementioned  carpet  plots  for  were  compared  to  another  figure 
in  the  same  Reference  (reproduced  here  as  Figure  33)  which  shows  the  lift 


loss  as  a  function  of  V^/V^.  This  question  is  discussed  at  the  end  of  Ap¬ 
pendix  A. 
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Effect  of  the  Stabilizer 

The  movable  horizontal  stabilizer  has  considerable  effects  on  moment 

and  lift.  Its  effect  on  drag  is  small  and  was  neglected. 

Again,  inconsistencies  in  the  data  of  Reference  5  were  noted.  Two  sets 

of  data  presented  there  do  not  agree .  The  aforementioned  carpet  plots  for 

C, ,  Cn,  and  C  give  values  at  two  positions  of  the  horizontal  stabilizer. 

L  D  in 

Therefore,  assuming  linearity,  C,  and  C  can  be  calculated.  Graphs  of 

X  \ 

C.  (V./V  ,  0.)  and  C  (a,  V./V  .  9.)  (Figures  37  to  42)  were  also  pre- 
L*  J  °°  J  m*  J  ®  J  r 

°H  H 

sented  in  the  same  reference.  When  the  derivatives  obtained  from  the  car¬ 


pet  plots  were  superimposed  upon  the  former  data  (at  8°  angle  of  attack. 
Figures  43  and  44) ,  the  carpet  plot  data  were  found  to  be  discontinuous  and 
inconsistent.  The  former  data  (Figures  37  through  42)  were  used  in  the 
HARRIER  program. 

However,  the  C,  and  C  data  given  do  not  include  data  for  V./V  =  0 
Lp  nu  j  00 

SH  SH 

or  0 .  =0.  Examination  of  Figure  44  reveals  that  C  is  close  to  linear 

3  SH 

with  both  V./V  and  9.,  and  linear  extrapolation  to  0.  =>  0  shows  that  C 
J  •  J  J  mS 

H 

at  any  V./V  converges  to  approximately  the  same  point  (-.0208).  Because  of 
J 

these  facts,  it  was  decided  that  linear  extrapolation  would  be  used  for  val¬ 
ues  of  0.  less  than  30°  and  for  values  of  V./V  less  than  5.4. 

J  J  00 

Examination  of  Figure  37  shows  that  Cr  (0.)  is  not  far  from  linear  for 

L5  j 
H 

Q.  between  30“  and  60°;  linear  extrapolation  for  values  of  9^.  less  than  30° 
was  used.  Although  the  change  in  CT  as  a  function  of  V./V  seems  to  be  a 

function  of  V./V  squared,  for  simplicity  linear  extrapolation  was  again 
J  00 

used  for  values  of  V./V  less  than  5.4. 

J  ® 
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The  data  in  Figures  37  through  42  were  used  for  Ct  and  C  in  the 

nig 

HARRIER  program.  H  H 

Inlet  Effects  Due  to  Ram  Drag 

The  engine  inlets  on  the  Harrier  are  considerably  forward  of  the  eg,  so 
ram  drag  at  the  inlet  affects  not  only  lift  and  drag  but  has  a  considerable 
effect  on  pitching  moment.  Ram  drag  as  a  function  of  engine  RPM  and  Mach 
number  was  taken  from  Reference  7  and  is  presented  in  Figure  45. 

The  change  in  pitching  moment  coefficient  due  to  ram  drag  was  computed 
in  Reference  S.  Referring  to  Figure  46, 

D 

C  »  ---?m  [CIO. 41  sin  (a  -  7.5°)  +  8.10  sin(7.5°  +  1.5°)] 

minlet  qSc  - - — — '  v ^ ' 

inlet  face  internal 

The  first  term  inside  the  brackets  represents  turning  the  flow  from 
freestream  to  the  inlet  face.  The  second  term  represents  turning  the  flow 
inside  the  inlet  (a  moment  arm,  in  feet ,  times  the  sine  of  the  angle  the 
flow  was  turned) . 

The  change  in  lift  coefficient  due  to  ram  drag  at  the  inlet  can  be  com¬ 
puted  in  a  similar  manner, 

_  q  sin(a  +  1.5a)  _  ^inlet 

Hnlet  Tam  qS  qS 

where  1.5°  is  the  additional  angle,  past  the  fuselage  reference  line,  by  which 
the  airflow  is  turned  to  the  engine  face.  The  above  equation  was  inserted 
into  the  HARRIER  program  in  the  form  of  Figures  47  and  48.  A  cross  check 
between  the  data  in  References  5  and  7  confirmed  the  inlet  effect  data. 
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Speedbrake  Effects 

As  given  by  Reference  5,  the  speedbrake  effects  are  functions  of  a  only. 
The  data  given  are  for  a  speedbrake  deflection  of  30°  and  all  analyses  were 
done  with  the  speedbrake  simulated  in  that  position.  Actual  speedbrake  pos¬ 
ition  with  undercarriage  down  is  fixed  at  25°. 

Figures  49  through  51  present  and 

speedbrake  speedbrake 

C 

“speedbrake 

Reaction  Control  System  (RCS) 

Bleed  air  from  the  high  pressure  stage  of  the  compressor  of  the 
Pegasus  11  engine  is  ducted  to  the  nose,  tail  and  wingtips  of  the  Harrier 
to  provide  moments  about  all  three  axes  at  low  speeds.  The  bleed  air  con¬ 
trol  is  phased  in  during  main  nozzle  deflection.  For  9^.  _>  20°  full  RCS 

moments  are  available.  For  0.  =  10°,  50%  of  maximum  moment  is  available. 

J 

The  RCS  is  not  operational  at  0^  =0°, 

The  moment  available  about  any  axis  (under  standard  day,  sea  level  con¬ 
ditions)  is  primarily  a  function  of  engine  speed  (RPM) ,  control  surface 
position,  and  total  demand  on  the  RCS.  The  maximum  amount  of  pitching  mo¬ 
ment  produced  by  the  pitch  control  thrusters  is  reduced  when  roll  and/or  yaw 
RCS  thrust  is  also  being  used.  Pitch  reaction  control  is  controlled  by  the 
longitudinal  stick,  roll  is  controlled  by  the  lateral  stick,  and  yaw  is  con¬ 
trolled  by  the  rudder  pedals.  Although  roll  and  yaw  reaction  controls  are 
not  used  in  this  study,  they  could  be  added  to  the  model. 

Pitch  control  comes  from  fore  and  aft  thrusters,  exhausting  downward, 
providing  nose-up  and  nose-down  moments,  respectively.  Figure  53  is  a  func¬ 
tional  block  diagram  showing  how  the  reaction  control  system  was  modelled. 
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Duct  pressure  ratios  (DPR  1  and  DPR  2)  are  computed  as  functions  of  engine 
RPM  and  of  total  reaction  control  system  nozzle  area  (TRCSA)  and  are  shown 
in  Figures  52  and  54.  Forward  and  aft  RCS  nozzle  areas (FRCSP  and  ARCSP)  are 
functions  of  the  stabilizer  position  (Figure  55).  The  percentage  of  thrust 
from  each  nozzle  (FRCST  and  ARCST)  is  then  a  function  of  nozzle  area  (FRCSA 
and  ARCSA,  Figure  56)  times  the  duct  pressure  ratios.  This  percentage  is 
then  multiplied  by  the  maximum  thrust  available  if  the  nozzle  were  full  open 
(FRCSMT  and  ARCSWT)  to  get  actual  thrust  at  the  nozzle.  This  thrust  is  fur¬ 
ther  multiplied  by  a  reaction  control  effectiveness  (DLRCV,  Fig.  57)  which  is 
a  function  of  velocity.  The  RCS  effectiveness  models  the  freestream  effects 
on  the  nozzle  exhaust.  The  freestream  tends  to  change  the  effective  angle 
of  the  nozzle  exhaust;  the  faster  the  aircraft  is  traveling  the  less  effec¬ 
tive  are  the  reaction  controls.  Figure  57  shows  RCS  effectiveness  (DLRCV) 
as  a  function  of  V  . 

oo 

Engine  Model 

The  Rolls  Royce  Pegasus  11  turbofan  engine  is  modelled  as  a  function  of 
static  gross  thrust  versus  revolutions  per  minute  (RPM).  The  engine  has  six 
distinct  thrust  ratings.  These  thrust  ratings  arise  from  temporarily  allow¬ 
ing  jet  pipe  and  turbine  over-temperatures ,  overspeeding  of  the  engine,  and 
injecting  water  into  the  engine  to  cool  the  combustion  chamber  and  turbine. 

At  each  increased  thrust  rating  there  is  a  maximum  operating  time  allowed. 
After  this  allowance  is  expended  the  engine  must  be  run  for  five  minutes  at 
or  below  the  maximum  continuous  rating  before  again  exceeding  it.  Each  en¬ 
gine  also  has  a  specified  engine  life  and  an  engine  life  counter  which  keeps 
track  of  time  and  temperature  when  operating  above  630°C  jet  pipe  temperature. 
Under  standard  day  conditions  this  is  equal  to  about  92%  RPM  or  a  T/W  of  1.01 
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(at  14000  lbs).  Therefore,  thrust  ahove  14000  lhs.  must  be  avoided  whenever 
possible  in  order  to  save  engine  life.  Figure  58  shows  gross  thrust  versus 
RPM  and  the  six  thrust  ratings.  Table  1  shows  the  equivalent  T/W  ratios  for 
a  14000  lb.  aircraft. 

Gross  thrust  increases  as  a  function  of  freestream  velocity  due  to  the 
higher  pressures  at  the  compressor  face.  This  gross  thrust  increase  (Fig¬ 
ure  60)  is  added  to  the  aforementioned  thrust  vs.  RPM  curve.  Knowing  the 
thrust  required  to  trim  the  aircraft,  the  RPM  can  then  be  found.  Given  RPM, 
Figure  59  shows  the  thrust  distribution  between  the  forward  and  aft  main 
nozzles  and  so  the  pitching  moment  as  a  function  of  thrust  is  easily  ob¬ 
tained. 

All  data  are  for  a  standard  50°F  day  at  sea  level.  The  effects  of  tem¬ 
perature,  humidity,  and  altitude  were  neglected.  The  effect  of  bleeding  RCS 
air  on  the  gross  thrust  was  also  neglected. 
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CHAPTER  III 


RESULTS  AND  DISCUSSION 


First  in  this  chapter  the  interpretation  of  trim-maps  is  discussed 
in  general.  This  is  followed  by  a  discussion  of  those  transition  con¬ 
straints  which  are  most  significant  for  the  AV-8A  Harrier.  General  con¬ 
clusions  concerning  significant  characteristics  of  the  AV-8A  transition 
are  then  presented.  Following,  there  is  an  example  using  the  trim-maps  to 
study  the  transition  profile  presently  recommended  for  use  during  shipboard 
landings.  The  chapter  is  followed  by  a  summary  of  the  most  significant 
conclusions . 

Since  landing  transitions  with  the  Harrier  are  more  difficult  than 
takeoff  transitions,  only  the  former  are  discussed  in  this  report.  However, 
trim-maps  for  a  climbing,  accelerating  transition  are  also  presented. 

Trim-maps  of  the  Harrier  in  transition  are  shown  in  Figures  65 
through  106*;  these  were  calculated  for  a  freestream  velocity  range  of  250 
ft/sec  to  0  ft/sec.  The  configuration  studied  was  14000  lbs  gross  weight 
and  7.8%  mac  eg.  Flaps  (<5p  =  50°)  and  undercarriage  were  down  and  speed- 
brake  was  out  at  30°  (normal  speedbrake  deflection  on  approach  is  actually 
25°).  The  HARRIER  also  can  handle  flaps-up  flight  with  the  under¬ 
carriage  up  and  the  speedbrake  still  extended,  with  any  weight  and  eg 
position.  Table  2  summarizes  the  flight  conditions  mapped  and  shows  where 
the  results  may  be  found. 

The  trim-map  formats  were  developed  by  T.  A.  Dukes,  G.  J.  Bom  and  T.  0. 

Williams  under  NASA  Grant  No.  NGR-31-001-277. 


Each  set  of  trim-maps  shows  families  of  curves  of  a(ALFA),  5„(DH), 

n 

T/W (THRUST) ,  L/W,  D/W,  9  (THETA  F)  and  V./V  as  functions  of  V  and  0 . 

t  3  0  00  3 

for  a  given  flight  path  angle  and  acceleration  along  the  flight  path. 

Each  plotted  point  represents  the  "trim"  value  of  the  variable  at  the 
given  6^  and  V^.  On  the  trim-maps  a  0  indicates  0  =  0;  a  "1"  indicates 
10  degrees;  a  "2"  20  degrees  of  9^  and  so  on.  An*  indicates  0^  =  100°. 

For  coincident  points  there  is  an  A  plotted  when  two  points  coincide;  a 
B  represents  three  coinciding  points,  etc.  The  data  points  of  all  the 
variables  are  truncated  whenever  a  exceeds  20°  or  T/W  exceeds  l.S. 

Creating  a  Flight  Profile 

It  must  be  emphasized  that  all  the  data  points  are  for  quasi-steady 
equilibrium  (trimmed)  flight  and  the  dynamics  involved  in  moving  from  one 
trim  point  to  another  are  not  considered  here.  However,  it  shall  be  seen 
that,  in  general,  moving  from  one  trim  point  to  another  is  a  very  straight¬ 
forward  process.  The  transition  profile  is  defined  as  all  the  parameters 
and  variables  which  represent  the  flight  path,  aircraft  configuration,  and 
flight  condition  at  which  the  aircraft  is  flying. 

Creating  a  transition  profile  means  piecing  together  the  trim  flight 
conditions  as  defined  in  the  trim-maps  in  order  to  create  a  desired  flight 
path  and  to  indicate  a  pilot  technique.  A  profile  may  be  created  by  de¬ 
fining  the  flight  path  to  be  taken  and  then  choosing  the  pilot  technique 
required  to  fly  that  path  or  by  defining  the  method  of  flying  and  working 
backward  to  determine  the  resultant  flight  path. 


The  desired  profile  is  constrained  by  many  factors,  some  of  which 


are  shown  directly  on  the  trim-maps.  For  example,  for  the  Harrier  the 
stabilizer  position  (<5^)  is  critical  because  there  is  only  marginal  nose- 
down  pitching  moment  available.  On  the  <S„  trim-maps  (Figures  66,  73,  80, 

n 

87,  94,  and  101)  it  can  be  seen  that  the  stabilizer  limit  of  +11.25 
degrees  is  exceeded  for  several  speeds  and  nozzle  positions.  This  means 
that  there  is  insufficient  nose^down  trim  moment  available  at  these  flight 
conditions,  not  even  considering  the  need  for  a  control  margin.  Similarly, 
angle  of  attack  (a)  and  thrust  (T)  limits  are  also  constraints  on  the 
transition  profile.  Obviously  the  maximum  gross  thrust  available  is  a 
limiting  factor  and  the  maximum  angle  of  attack  should  be  kept  below  buffet 
(16-18°).  Thrust  impact  on  engine  life  and  fuel  consumption  are  also  con¬ 
siderations.  Fuel  consumption  depends  on  gross  thrust,  which  in  turn 
depends  on  the  deceleration  profile  to  hover  (Reference  8).  Therefore, 
accurate  flight  path  control  is  not  only  important  for  safety  reasons  but 
is  also  important  for  fuel  considerations.  Pilot  visibility  makes  it 
desirable  that  the  aircraft  maintain  a  relatively  low  pitch  attitude  so- 
that  the  pilot  has  a  good  view  over  the  nose. 

In  a  landing  transition,  deceleration  and  descent  trim-maps  are 
generally  of  the  greatest  interest.  It  should  be  noted  that  a  deceleration 
of  0.1  g  with  the  flight  path  angle  (Y)  equal  to  zero  results  in  essentially 
the  same  trim-maps  as  zero  deceleration  and  Y  =  -0.1  radian  because  the 
component  of  weight  along  the  flight  path  in  a  descent  of  y  =  -0.1  rad  is 
equal  to  0.1  times  the  aircraft  weight.  This  equivalence  is  demonstrated 
by  Figures  79  through  85  and  Figures  86  through  92  which  show  the  trim-maps 
for  y  =  -0.1  rad,  uy  =  0.0  and  Y  =  0.0,  uy  =  -O.lg  respeetively .  There¬ 
fore,  for  small  angles  the  sum  of  the  flight  path  angle  in  radians  and  the 
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acceleration  in  g's  is  an  important  characteristic  of  the  trim-maps.  For 
example,  a  set  of  y  =  -0.1  rad.,  uy  =  -0.2  g  trim-maps  is  nearly  identical 
to  the  Y  =  -0.2  rad.,  uy  =  -0.1  g  (sum  =  -0.3)  trim-maps,  therefore  only 
one  of  the  two  sets  needs  to  be  computed. 

Control  Moment  Available 

The  critical  speeds  in  the  transition  of  the  AV-8A  Harrier  are  between 
approximately  225  ft/sec  and  12S  ft/sec  because  of  the  marginal  nose-down 

control  moment  available  in  this  range.  At  the  higher  speeds  the  control 

moment  is  provided  mostly  by  the  aerodynamic  force  acting  on  the  horizontal 
tail.  At  low  speeds  the  pitching  moment  is  produced  by  the  reaction  control 
system  which  looses  effectiveness  at  higher  speeds  (Figure  57)  because  the 
freestream  deflects  the  RCS  exhaust  and  changes  the  effective  exhaust  angle. 
In  the  225  to  125  ft/sec  speed  range  the  aerodynamic  pitching  moment  due 
to  stabilizer  deflection  is  decreasing  while  the  RCS  pitching  moment  is  not 
yet  strong  enough  to  make  up  for  the  loss.  Figures  61  and  62  demonstrate 

this  problem.  At  250  ft/sec  the  stabilizer  produces  almost  2400  ft-lb 

per  degree  of  5^  deflection.  With  the  trailing  edge  full  down  (11.25°) 
this  results  in  a  pitch-down  moment  of  almost  12600  ft-lb.  This  decreases 
to  less  than  4725  ft-lb  at  150  ft/sec,  a  loss  of  about  7875  ft-lb  of 
aerodynamic  pitching  moment.  Figure  62  shows  that  the  pitching  moment 
available  from  the  reaction  control  system  increases  from  about  1700  ft-lb 
at  250  ft/sec  to  4400  ft-lb  at  150  ft/sec.  an  increase  of  2700  ft-lb. 
Therefore,  slowing  down  from  250  to  150  ft/sec,  there  is  a  net  loss  of 
over  5000  ft-lb  in  the  available  pitching  moment. 

Two  critical  parameters  affecting  the  required  trim  moment  in  this 
speed  range  are  the  speedbrake  deflection  and  the  effects  of  thrust. 
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Speedbrake  Effects 

Without  the  speedbrake  extended  the  aircraft  exhibits  an  undesirable 
oscillatory  roll  mode  at  slow  speeds.  The  speedbrake  is  extended  on 
approach  not  for  its  drag  but  because  it  alters  the  airflow  around  the 
ventral  fin  (Reference  9) . 

As  mentioned  before,  the  limited  nose-down  pitching  moment  available 
is  a  quite  stringent  constraint  on  how  the  transition  may  be  flown  with 
the  AV-8A.  A  speedbrake  deflection  of  25°  (the  data  used  in  this  report 
was  for  30°  deflection)  causes  an  unfavorable  nose-up  aerodynamic  pitching 
moment  equivalent  to  about  3°  stabilizer  deflection  at  a  =  4°  and  1° 
stabilizer  deflection  at  oi  =  8°  (Figure  51).  With  the  speedbrake  fully 
extended  (66°)  the  ACm  due  to  speedbrake  is  nose  down  approximately  -0.02, 
(Reference  S,  flaps  §  U/C  up,  0^  =  0)  and  equivalent  to  about  one  degree 
stabilizer  deflection.  Therefore  it  would  be  desirable  to  use  full  speed- 
brake  deflection  during  transition  if  it  could  be  established  that  this 
would  stabilize  the  oscillatory  roll  as  well  as  the  25°  deflection.  As 
the  aircraft  was  decelerated  to  below  the  critical  speed  range  the  speed- 
brake  would  then  be  retracted  back  to  25°  to  keep  it  from  hitting  the 

i 

ground  at  touchdown. 

Thrust  Effects 

The  pitching  moment  varies  significantly  as  a  function  of  the  main 
nozzle  thrust,  T.  At  low  thrust  the  bypass  airflow  (fan  exhaust)  produces 
more  thrus^t  than  the  turbine  exhaust.  Since  the  fan  exhaust  flows  out 
the  front  nozzles  and  the  turbine  exhaust  out  the  aft  nozzles,  a  significant 
nose  up  moment  is  created.  Figure  59  shows  that  the  ratio  of  fan/turbine 
thrust  decreases  as  the  total  thrust  (or  RPM)  increases,  and  as  a  result 
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the  nose-up  moment  changes  to  a  nose-down  moment.  The  variation  of  the 
pitching  moment  with  thrust  is  shown  in  Figure  63.  Because  of  the  marginal 
nose-down  moment  available  it  would  be  advantageous  to  run  the  engine  at  a 
high  thrust  to  get  the  benefit  of  the  nose-down  moment  due  to  thrust. 
Running  the  engine  at  higher  RPM  also  increases  the  thrust  of  the  RCS, 
enabling  the  aft  pitch  thruster  to  produce  more  nose-down  moment  than  at 
low  thrust  settings. 

While  running  the  engine  at  higher  thrust  levels  improves  the  control 
problem  it  conflicts  with  one  objective  of  landing;  that  is,  to  descend  and 
decelerate  using  as  little  fuel  as.  possible.  Using  thrust  greater  than 
about  92%  also  decreases  engine  life.  This  should  be  avoided  in  landing 
transitions  as  much  as  possible  in  order  to  conserve  engine  life  for  take¬ 
off  and  combat.  Nevertheless,  a  trade-off  between  fuel  economy  and  an 
improved  control  margin  is  possible. 

Miminum  Fuel  Consumption 

The  thrust  trim-maps  indicate  implicitly  the  trim  points  for  minimum 
fuel  consumption  flight.  They  are  simply  the  8j  at  each  airspeed  requiring 
the  least  thrust.  The  minimum  thrust  points  in  the  trim-maps  (Figs.  67, 

74,  81,  88,  95,  102)  are  determined  by  the  a  =  20°  truncation.  Allowing 
for  an  a  margin  to  stall,  a  much  lower  value  of  a  must  be  cross-plotted 
on  the  thrust  curves.  As  an  example,  in  Figure  95  the  ot  =  12°  values  for 
various  0^  are  indicated  by  X-es;  in  the  decelerating-descending  flight 
of  that  Figure,  the  X-es  denote  the  truncation  points  of  the  thrust  curves 
for  a  *  12°.  A  minimum  thrust  program  with  this  a  constraint  amounts  to 
flying  at  this  angle  of  attack  essentially  throughout  the  transition.  In 
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order  to  establish  the  control  deflections  needed  for  such  a  transition, 
X-es  corresponding  to  a  =  12°  are  also  shown  in  the  6  curves  in  Figure  94. 
The  minimum  thrust  point  at  9^  =  80°  clearly  would  not  leave  an  adequate 
control  margin.  This  graph  indicates  that  below  about  175  ft/sec  the 
minimum  thrust  policy  should  be  abandoned  in  favor  of  an  improved  control 
margin.  This  illustrates  how  the  trim-maps  can  be  used  in  constructing 
transitions,  combining  requirements  and  constraints. 

Sample  Transition  Profile 

In  this  section  the  trim-maps  are  used  to  study  a  particular  transi¬ 
tion.  In  this  case  the  flight  program  is  known  and  the  trim-maps  will 
be  used  to  illustrate  how  the  aircraft  is  flown  and  the  pilot  technique 
involved. 

The  sample  transition  profile  is  the  one  presently  being  adopted  for 
use  by  the  United  States  Marine  Corps  for  shipboard  operations.  This 
profile  is  shown  in  Figure  64.  Beginning  at  12  nautical  miles  (nm)  the 
aircraft  is  at  1200  ft  ,  340  ft/sec,  flaps  and  undercarriage  down,  and 
nozzles  aft.  The  aircraft  is  then  slowed  to  approximately  290  ft/sec  and 
the  nozzles  rotated  to  20°  by  the  time  it  is  8  nm  from  the  ship.  Then 
from  8  to  4  nm  the  pilot  flies  a  constant-speed  descent  to  about  450  ft 
altitude  and  levels  out.  At  4  nm  and  270  ft/sec  the  pilot  puts  the  nozzles 
down  to  40°  and  further  decelerates  to  about  240  ft/sec  by  the  time  he  is 
1  nm  from  the  ship  (still  at  450  ft).  At  this  point  the  final  descent 
is  begun.  The  nozzles  are  lowered  to  81°  and  the  pilot  simultaneously 
decelerates  and  descends  to  a  hover  above  the  ship. 
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The  final  mile  of  the  transition  is  by  far  the  most  challenging  so 


it  will  be  traced  on  the  trim-maps. 

At  1  nm: 

=  240  ft/sec 

0.  =40° 

J 

h  =  450  ft 

On  Figures  72  through  79  point  "a"  represents  the  aircraft  flight 
conditions  for  straight  and  level,  unaccelerated,  flight  at  240  ft/sec  and 
9^  =  40°.  From  these  trim-maps  it  can  be  seen  that 

at  point  a,  =  240  ft/sec,  uy  =  0.0  g,  y  =  0.0  rad,  0  =  40°: 

a  =  9°;  SH  =  3°;  0  =  9°; 

T/W  =  0.48;  L/W  =  0.65;  D/W  =  0.205;  V./V  =  0.28 

j  00 

The  angle  of  attack  found  on  the  trim-map  agrees  closely  with  the  actual 
approach  angle  of  attack  in  that  the  Marine  Corps  pilots  attempt  to  keep 
a  constant  at  8°. 

At  this  point  the  nozzles  are  rotated  down  to  80  degrees  and  the 
pilot  begins  a  decelerating  descent.  If  the  pilot  held  the  angle  of 
attack  constant  while  he  rotated  the  nozzles  the  aircraft  would  begin  to 
decelerate  at  6.4  ft/sec^  (0,2  g)  and  descend  at  y  =  -5.7°  (-0.1  rad) 

(as  shown  by  point  b  in  figures  100  through  106).  Summarizing  these 
trim-maps : 
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at  point  b,  =  240  ft/sec,  uy  =  -0.2g,  y  =  -0.1  rad,  0  =  80° 

a  =  9°,  6h  =  4.75° 

T/W  =  0.47,  L/W  =  0.S5,  D/W  =  0.195 

The  trim-maps  verify  that  this  is  a  good  place  to  rotate  the  nozzles 
from  40°  to  80°  because  a  and  thrust  remain  essentially  unchanged.  There¬ 
fore,  besides  rotating  the  nozzles, the  pilot  only  needs  to  move  the  stick 
forward  enough  to  move  the  elevator  1.75°  to  change  from  one  quasi¬ 
equilibrium  trim  point  to  another. 

Examining  further  Figures  100  through  106  it  can  be  seen  that 
maintaining  the  -0.2  g  deceleration  and  -0.1  rad  flight  path  angle  requires 
increasing  a  and  5U,  and  at  160  ft/sec  a  would  be  16°  and  S„  would  be 
10.75°,  both  dangerously  high.  Angle  of  attack  would  be  approaching  stall 
and  6^  would  be  approaching  the  forward  stick  stop  (11.25°).  Examining 
the  trim-maps  for  unaccelerated  straight  and  level  flight  (Figures  72 
through  78)  and  for  decelerating  descent  at  uy  =  -0.1  g,  y  =  -0.1  rad 
(Figures  93  -  99),  it  can  be  seen  that  the  flatter  the  flight  path  and 
the  less  the  deceleration,  the  less  the  stabilizer  required.  This  is 
caused  mainly  by  the  fact  that  the  trim  condition  for  lower  deceleration 
requires  more  thrust  (see  Figs.  95  and  102),  which  decreases  the  nose  up 
moment  due  to  thrust. 

As  a  result,  assuming  the  maximum  a  allowed  to  be  about  12.5°,  the 
uy  =  -.2  g,  y  =  -0.1  rad  flight  may  only  be  maintained  until  200  ft/sec 
is  reached  (point  c)  and  the  deceleration/descent  must  then  be  decreased. 
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Referring  to  Figures  100  -  106: 


at  point  c,  Vm  =  200  ft/sec,  uy  =  -0.2  g,  Y  =  -0.1  rad,  0^  =  80° 

a  =  12.5°;  6U  =  6.25°;  0=7° 


T/W  =  0.55;  L/W  =  0.46;  V./V  =  0.37;  D/W  =  0.175 

J  00 


To  fly  from  the  conditions  at  point  b  to  point  c  the  pilot  would 
simply  continue  to  move  the  stick  slowly  forward  to  6.25°  stabilizer  angle, 
allowing  the  angle  of  attack  to  increase  to  12.5°.  During  this  time  thrust 
must  be  increased  by  a  small  amount  (0.08  T/W). 

At  this  point  the  angle  of  attack  must  either  be  held  constant  or  be 
reduced.  This  means  the  nozzle  deflection  must  either  be  increased  or  the 


deceleration/descent  reduced.  Since  the  Marine  Corps  pilots  do  not  move 
the  nozzles  past  82°  (hover  stop)  it  is  therefore  evident  that  the  de¬ 
celeration/descent  angle  is  reduced.  In  addition,  the  pilots  are  constantly 
trying  to  keep  the  angle  of  attack  near  8°  and  as  hover  is  approached  they 
decrease  the  deceleration  while  flattening  the  flight  path. 

This  means  the  pilot  is  trying  to  fly  smoothly  between  the  a  =  8° 
trim  points  on  the  uy  =  -0.1  g,  y  =  -0.1  rad  trim-maps  (Figures  93-99)  and 
the  uy  =  0.0  g  and  y  =  0.0  rad  trim-maps  (Figures  72-78).  From  200  ft/sec 
to  100  ft/sec,  the  pilot  should  take  special  care  to  keep  the  angle  of 
attack  at  8°  or  above  in  order  that  the  pitch-up  due  to  the  speedbrake 
(Figure  51)  is  avoided.  This  speed  range  can  be  traced  on  the  trim-maps 
in  the  following  way.  Comparing  Figures  93  and  72,  note  that  points  d  and 


e  represent  the  8  =  80°,  a  =  8°  trim  points  on  the  a  trim-maps  for 

•  • 

uy  ■  -0.1  g,  Y  *  -0.1  rad,  and  for  uy  =  0.0  g,  Y  =  0.0  rad,  respectively. 


Point  d  shows  trimmed  decelerating,  descending  flight  at  185  ft/sec  at  an 


angle  of  attack  of  8°.  As  the  deceleration  and  descent  angle  are  decreased, 
the  airspeed  for  trimmed  flight  at  a  =  8°  with  the  same  0^  is  decreased  to 
50  ft/sec  when  steady  level  flight  is  reached  (point  e) .  The  airpseed  for 
trimmed  flight  at  8°  for  uy  =  0.05  g,  Y  =  -0.05  rad  can  be  estimated  to  be 
about  118  ft/sec  [(185  +  50) /2], 

Using  this  information,  the  pilot  technique  and  flight  conditions 
for  flying  this  path  can  be  deduced.  From,  the  conditions  given  at  point  c, 
the  stick  will  have  to  be  moved  to  almost  full  forward  at  about  200  ft/sec 
as  the  amount  of  deceleration  is  decreased  to  -0.1  g  (point  d,  Fig.  94). 

As  the  deceleration  and  flight  path  angle  are  further  decreased,  the  stick 
will  be  very  slowly  moved  aft  until  about  150  ft/sec  is  reached  when  it 
sould  be  moved  more  rapidly  aft  in  the  150  ft/sec  to  75  ft/sec  speed  range. 
In  the  75  ft/sec  to  50  ft/sec  speed  range,  the  stick  would  continue  to  be 
moved  slightly  aft  to  about  6  =  3.25°.  During  this  time,  thrust  would 

be  initially  rapidly  increased  to  about  T/W  =  0.8  and  then  more  slowly 
increased  reaching  T/W  «  1  at  about  100  ft/sec. 

After  level  unaccelerated  flight  is  reached  at  50  ft/sec  (point  e) 
simply  raising  the  angle  of  attack  to  10°  will  slowly  bring  the  aircraft 
to  hover  (point  f  on  Figures  72-78) . 

This  example  illustrates  how  the  trim-maps  can  be  used  for  a  detailed 
analysis  of  a  prescribed  transition.  In  the  AV-8A  the  control  power  limita¬ 
tion  and  the  simplest  possible  piloting  technique  are  overriding  factors. 

In  future  VTOL  aircraft  with  no  severe  control  limitation  and  with  various 
degrees  of  automatic  flight  control,  the  trim-maps  can  be  used  to  explore 
less  constrained  and  potentially  more  efficient  transition  profiles. 


CONCLUSIONS 


A  mathematical  computer  model,  "HARRIER",  of  the  AV-8A  Harrier 
aircraft  has  been  developed  utilizing  the  concept  of  quasi-stationary  flight 
in  order  to  determine  quasi -equilibrium  values  of  flight  variables  and 
controls  during  accelerating-decelerating  flights.  Using  available  aero¬ 
dynamic  data  via  look-up  schemes,  "trim-maps”  are  computed  which  show  the 
values  of  the  longitudinal  flight  and  control  variables  as  functions  of  air¬ 
speed  and  thrust  angle.  These  maps  have  been  shown  to  be  useful  tools  in 
piecing  together  transition  profiles,  taking  into  consideration  various 
requirements  and  constraints. 

The  AV-8A  is  most  difficult  to  control  longitudinally  in  the  speed 
range  of  approximately  22S  to  12S  ft/sec.  Decelerating  through  this  critical 
speed  range  the  aerodynamic  pitching  moment  from  the  stabilizer  is  decreas¬ 
ing  as  a  function  of  dynamic  pressure  and  the  reaction  control  system  has 
not  yet  become  effective  enough  to  make  up  for  the  loss.  Considering  that 
the  transition  is  constrained  by  the  available  nose-down  control  moment,  the 
moments  caused  by  the  speedbrake  and  the  thrust  become  significant.  There¬ 
fore,  the  transition  should  be  flown  with  two  major  constraints  in  this 
speed  range.  First  the  angle  of  attack  should  be  kept  at  or  above  8°  so 
that  the  nose-up  pitching  moment  due  to  the  speedbrake  is  not  excessive. 
Second,  deceleration  through  the  critical  speed  range  should  be  performed 
gradually  to  keep  thrust  up  near  the  thrust  required  for  level  flight  and 
therefore  minimizing  the  nose  up  pitching  moment  due  to  thrust.  This  second 
constraint  should  be  tempered,  especially  at  higher  landing  weights,  with 
the  consideration  of  engine  life. 


Extending  the  speedbrake  to  the  full-out  position  on  approach  would 
be  beneficial  in  that  it  would  change  the  pitch-up  effect  of  the  speed- 
brake  to  pitch-down,  thereby  alleviating  the  problem  of  marginal  nose  down 
moment  available.  After  decelerating  below  100  ft/sec  the  speedbrake  should 
then  be  retracted  to  25°  to  avoid  hitting  the  ground  at  touchdown. 

The  final  decelerating-descending  phase  of  the  present  USMC  approach 
profile  has  been  analyzed  in  detail  with  the  aid  of  the  trim-maps.  The 
technique  of  flying  the  transition  as  suggested  by  the  trim-maps  has  been 
found  in  good  agreement  with  the  techniques  actually  used  under  the  given 
constraints.  This  established  confidence  in  both  the  mathematical  model 
and  in  the  trim-maps  as  tools  for  studying  alternate  transition  profiles 
which  might  improve  pilot  work  load  and/or  fuel  efficiency. 
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TABLE  1 


Thrust  to  Weight  Ratios  and  Time  Limitations  of 
The  Installed  Pegasus  11  Engine* 


THRUST 

% 

TIME 

•GROSS 

T/W 

RATING 

nf 

LIMITATION 

THRUST 
**  165~ 

8  14000  lbs 

Max.  continuous 

90 

unlimited 

12900 

0.92 

Maximum  thrust 

94.4 

15  min. 

15500 

1.10 

Normal  lift,  dry 

100.0 

2.5  min. 

17700 

1.26 

Short  lift,  dry 

102.7 

15  sec. 

18650 

1.33 

Normal  lift,  wet 

104.1 

1.5  min. 

19300 

1.38 

Short  lift,  wet 

105.8 

15  sec. 

19600 

1.40 

•Static,  Sea  Level  Standard  Day 


TABLE  2 

Summary  of  Trim-Map  Flight  Conditions 


Figures 

Y 

• 

uv 

Cg) 

"v 

Cg) 

v. 

ft/sec 

0. 

J 

Cdeg) 

5f 

Cdeg) 

65-71 

0.1 

0.1 

0 

250-0 

0-100 

50 

72-78 

0 

0 

0 

250-0 

0-100 

50 

79-85 

-0.1 

0 

0 

250-0 

0-100 

50 

86-92 

0 

-0.1 

0 

250-0 

0-100 

50 

93-99 

-0.1 

-0.1 

0 

250-0 

0-100 

50 

100-106 

-0.1 

Cn| 

• 

o 

I 

0 

250-0 

0-100 

50 
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L  3  0.90W 

D  »  0.375W 
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a  »  7.0° 
9.  -  50° 


T  *  0.57W 
L  ■  0.53W 
D  »  0.190W 


FIGURE  1  -  Two  Sets  of  Control  Positions  Resulting  in  the  Same  Trimmed 
Flight  Condition,  (Level  Unaccelerated  Flight  »  220  ft/s). 
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FIGURE  5  -  AV-8  C.  With  and  Without  Wing/Jet  FIGURE  6  -  AV-8A  C.  With  and  Without  Wing/Jet 
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FIGURE  23  -  AV-8A  CQ  With  Wing/Jet  Interaction,  FIGURE  24  -  AV-8A  CD  With  Wing/Jet  Interaction 
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FIGURE  34  -  CL 
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Flaps  Up,  Without  Wing/ Jet  Interaction 
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FIGURE  35  -  CD  Flaps  Up,  Without  Wing/Jet  Interaction 

wing-body 

(Ref.  5,  p.  160). 
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FIGURE  36  -  C  Flaps  Up,  Without  Wing/ Jet  Interaction 

rawing-body 


FIGURE  45  -  Ram  Drag  at  Various  Mach  Numbers. 
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FIGURE  46  -  Modelling  AC.  Due  to  Ram  Drag  (Ref.  5). 
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FIGURE  47  -  Change  in  Lift  Coefficient  Due  to  Ram  Drag  FIGURE  48  -  Change  in  Lift  Coefficient  Due  to  Ram  Drag 

(AC,  )  for  Various  V.  RPM  =  94.4%.  (AC  )  for  Various  V  ,  RPM  =  105.8%. 
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FIGURE  53  -  Functional  Block  Diagram  of  the  AV-8A  Reaction  Control  System 
and  Engine. 
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FIGURE  55  -  RCS  Nozzle  Opening  as  a  Function  of 
Stabilizer  Position. 


FIGURE  56  -  Forward  and  Aft  RCS  Thrust  Ratios  as  Functions  of  Their 
Respective  Nozzle  Openings. 
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FIGURE  61  -  Pitching  Moment  per  Degree  of  FIGURE  62  -  Maximum  Nose  Down  RCS  Control 

Stabilizer  Deflection.  Pitching  Moment. 


FIGURE  63  -  Variation  of  Pitching  Moment  with  Gross  Thrust  at  0^  «  90°. 


FIGURE  64  -  Shipboard  Approach  Profile  as  Used  by  the  United  States 
Marine  Corps. 
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FIGURE  65  -  Trlimap  of  a(9,,  V  );  with  6, 
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FIGURE  69  -  Trlnmap  of  D/W(9  ,  V^);  with  <5| 
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FIGURE  70  -  Triraiap  of  8(9,,  V  );  wit1!  6, 
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FIGURE  71  -  Trinmap  of  V./V  (9  ,  V  );  with  6 
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HA»BIER;»  =  0-300;  TH  SI  J=0- Qfl.  8;  ?LAP  =  <50 . 0 ;  \  P? &  = 


FIGURE  72  -  Trlnmap  of  a(8.,  V  );  with  «_  =  50°,  0  »  0.0  g,  yr  *  0.0  rad. 
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(PT/SEC)  35 ;HHSRIEP;¥= 0-300;  THETJ=0- 98.5 ; PU ?=50 . 0; A  PH fi  = 

25J.000  - - 1 - 0 - I---1 - 


FIGURE  73  -  Trlmmap  of  «H(0 . .  VJ;  with  6F  =  50°,  0y  »  0.0  g,  yQ  =  0.0  rad. 


(FT/SBC)  35;  HARBieR  ;¥*!>- 300;  TH  ET.J=  0-OR. 5;  FLA  P  =  5C.  0;  Ap  A>-  = 


FIGURE  74  -  Trlnwap  of  T/W(8.,  VJ;  with  «F  *  50*.  0y  *  0.0  g,  yc  *  0.0  rad. 


|FT/SEC>  JS;M»PmEF;»»0-300:  TIHtJ'O-'JR.S;  »L».P  =  5’) 


FIGURE  75  -  Trlmnap  of  L/W(9.,  VJ;  with  6, 


IPT/SBC)  35;HdPRtEfi;I®0*3O0;  TH El J= 0- P0 . 5 ; PL  A P= 


FIGURE  76  -  Triwmap  of  O/W(0j,  VJ;  with  6, 


JPT/SBC)  3*5;  HARRIER;  V=  0-300;  rHET.J=C-'tfl  .  5  ;  FLM‘  =  50 

250.030  . . - . I . « . I  7— 


FIGURE  77  -  Trlimap  of  e(0,.  VJ ;  with  6p  =  50°,  0y  *  0.0  g,  yc  »  0.0  rad. 


(Fl/SCC)  5«;t»BBItB;»-0-300;  THElJ=0-98. 5; PLAP=50 . 0; APAB=  0.0;UAB= 


FIGURE  79  -  Triimtap  of  a(9  ,  V  );  with  «F  =  50°,  i)  *  0.0  g,  yc  ■  -0.1  rad. 


(Fl/SECI  5«;t4RRI£R;»*0-J00;  TH*lJ=0-98.5i FLAP=50. 0; APAR=  0.0;G*H=-0 

iSw.OOO  - - 2A--A - 5 - 6—1 - 7 - |  — 


FIGURE  81  -  Trlntnap  of  T/W(e.,  V  );  with  6, 
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(PT/SiC)  5«; BARRI18;V=0-300;  TH£lJ=Q-98. 5 ; FLA P=50 .0 ;A PAH=  0.0;GAN=-0 


FIGURE  84  -  Trlmnap  of  e(e.,  V  );  with  6. 


FIGURE  85  -  THirmap  of  V./V  (0.,  V  ?;  wit h  5,  «  50°,  0  *  0.0  g,  yr  *  -0.1  rad 
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(FI/S£C)  65;U»Rai.fcR;V=0-JOO;  THKTJ  =  0- 98.  5;  e'UP'SO .  0 ;  APAR  =  -0. 


FIGURE  86  -  Trlmmap  of  a(8j,  V„) ;  with  6| 


(rr/icC)  55; B»«RIEH;»=0-300;  TH ETJ*0-98 .5 ; ?LA P=50 . 0 ;A PA R=-0.  1;GArt= 


FIGURE  87  -  Trlntnap  of  6H(0, ,  V^) ;  with  6^  =  50°,  ■  -0.1  g,  *  0.0  rad. 


(Fl/SAC)  55;BARHIEB;t*Q-JOO;  TU*1J=0- 98.5 ; FLAP=50 . 0; AP A R=-0 . 1 ; GA«= 


FIGURE  88  -  Trinwap  of  T/W{e.,  VJ ;  with  •  50* ,  0  *  -0.1  g,  y  =  0.0  rad. 


(ft/SBC)  55:l»8HIER;»*0-3C0;  TH 114=0 -98 .5; FL»P=54. C; AP *8=-0. I ; GAH= 


FIGURE  90  -  Triiwap  of  0/W(8.,  VJ;  with  =  50’,  0  =  -0.1  g,  y  =  0.0  rad. 


(PT/S»C|  5S;E»BFIfH;V«C-300;  THETJ=0-98.S ; PU P=SO . 0; A PAB=-0 . 1 ;GAN= 


FIGURE  91  -  Trlnwap  of  0(0,,  V  ) ;  with  6  ■  50',  i)  *  -0.1  g,  Y  =  0.0  rad. 


(FT/SEC)  51 ; HARRIED ;f *0-300;  THETJ=0- 90. 5; FL»P=50 . 0; APAR=-0. 1;GAN=-0 


t 


FIGURE  93  -  Trinmap  of  a(0,,  V  );  with 
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(f* I/SEC)  5 1 ;  HAPRI8F  ;*»0* 300;  THETJ=0- 9fl. 5 ; PLAP=50. 0; »P*B=-0. 1 ; G*«=-0 


FIGURE  94  -  Trlnwap  of  6U(9.,  V  );  with  6, 
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(FT/SEC)  5 1 ; HAFRI ER ; F*0- 300 :  TH ETJ=C- 98.5 ; PLAP=50 . 0; APA R*-0. »;GA!1=-0 

25J.000  ♦ . 6 - - 7 . | . . . 8|-- 

I  X  I  I  • 

1  3A  511  G  0  7  |  8  | 


FIGURE  95  -  Trlnmap  of  T/W(9.,  V  J ;  wi th  6p  *  50* ,  0y  >  -0.1  9.  T£  !  -0.1  rad 


T/SECJ  51 ; HARRIER ;»*0-300;  TH FTJ= 0-98 .5 ; PIAP=50. 0; AP AR=-0. 1 ; GAS= 


(PT/StC)  51;  HU  BRIER;  V=  0-300;  TH£TJ=0-‘)8.5;FLAP=50.0;APfcR  =  '0. 1  ;GAN 


51;HARRIBB;*=0-300:  THBT J=0-9H . 5 ; FLAP=50 . 0 ; A PAR=-0.  1;GAR= 
— . 3B506--7-I  ---0 . I - - 1 


FIGURE  99  -  Trlnrup  of  V./VJe.,  VJ;  with  6, 


FIGURE  100  -  Trinmap  of  a(8.,  V  J ;  with  6p  *  50',  0y  «  -0.2  g,  y  *  -0.1  rad. 


(FT/SBC)  4 7 ;H»RRI ER ;9*0- 300;  THETJ=0- 96.  5 ;ELAP*50. 0 ; AP*fi  =  -0 .2;GAB=-0 


FIGURE  101  -  Trlmmap  of  6^(9.,  v„);  with  6^  *  50',  0^  *  -0.2  g,  y c  *  -0.1  rad. 
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FIGURE  103  -  Trlnroap  of  l/W(6.,  VJ;  with  «F  *  50°,  0y  =  -0.2 


(FT/SEC)  *7; HARRIER ;f«0- JOO;  TH fcTJ«0-98. 5; PL»P=50. 0;AP AR *-0. 2 ;G*K=-0 


FIGURE  104  -  Trlmnap  of  D/W(6.,  V  );  with  5,  -  50*.  d  =  -0.2  g,  y  *  -0.1  rad. 


FIGURE  105  -  Trlnnap  of  8(8.,  VJ;  with 


(fT/SEC)  *7;  HA  RRiea;¥  =0-300;  TH  E1J  *0 -98 . 5;  FLAP*  50. 0;AP  A8*-0. 2  ;GAH=-0 

250.000  ♦ - 01-7 - 8-| . . . . I-- . . . |-- 


FIGURE  106  -  Trlmnap  of  V./V  (0.,  V  );  with  6C  -  50*.  0  *  -0.2  g,  y  *  -0.1  rad. 
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APPENDIX  A 


EQUATIONS  OF  MOTION 

In  this  section  the  equations  used  in  modeling  the  transition  of  the 
AV-8A  will  be  presented.  First  the  equations  will  be  derived  for  three 
dimensional  flight,  one  set  at  a  time.  Then  each  set  of  equations  will  be 
reduced  to  consider  only  flight  in  a  vertical  plane  with  zero  sideslip.  As 
stated  in  Chapter  2,  all  aircraft  body  rotational  velocities  and  accelerations 
are  assumed  to  be  zero  in  quasi-equilibrium  flight.  Therefore,  Pg  =  qg  = 

rB  =  =  =  rB  °* 


TRANSLATIONAL  ACCELERATIONS 

Now  the  equations  for  tig,  vg  and  Wg  will  be  derived.  It  is  first 
assumed  that  the  rate  of  fuel  consumption  (jJ  has  a  negligible  effect  on  the 
dynamics  of  the  aircraft.  However,  the  fuel  consumed  throughout  a  maneuver 
may  still  be  considered  by  adjusting  W  to  reflect  fuel  consumed.  With  the 

assumption  W  =  constant. 

•  W  • 

F  =  mV  =  -  V 

g 


or 


•  F 

V  =  a  =  g-  for  acceleration  in  g's. 
w 

In  aircraft  body  axes 

is  =  (^B  +  Ib)/w  +  *B  (in  g's) 

where 


(A-l) 


A 


translational  accelerations,  body  frame  (in  g's) 


£g  A  y g  aerodynamic  forces,  body  frame 


Ib  ^ 


engine  and  reaction  control  system  forces,  body  frame 


gD  A  g  local  unit  gravity  vector,  body  frame 

^  -  yB 


While  is  normally  defined  in  body  axes,  the  aerodynamic  forces 

— D 

are  usually  expressed  in  vind  axes  and  the  gravity  force  in  the  local 

inertial  axes  (flat,  non-rotating  earth).  Therefore,  these  forces  must 

be  transformed  to  body  axes. 

F_  =  hB  f 
— B  w  -w 


F  (cos  a  cos  6)  (-  cos  a  sin  6)  (-  sin  a)  F 

XB  w 


(sin  3) 


(cos  3) 


(sin  a  cos  3)  (-  sin  a  sin  3)  (cos  a) 
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h 


H 


B 

I  S.I 


V 

B 

\ 

S3 

S_ 

( cos  9  cos  0) 


(sin  0  sin  9  cos  0 
-  cos  0  sin  0) 

(cos  (f  sin  9  cos  0 
+  sin  0  sin  0  ) 


(cos  9  sin  0) 


(sin  <(i  sin  9  sin  0 
+  COS  0  cos  0) 

(cos  0  sin  0  sin  0 
-  sin  0  cos  0) 


expanding 


(-  sin  9) 
(sin  0  cos  6). 


(bos  0  cos  93 


0 


1 


for  the  acceleration  of  gravity  in  g's 

a°d  Zg  rausi  be  determined.  includes  all  of  the  aerodynamic 
forces  acting  on  the  aircraft  including  effects  of  the  thrust  on  the  wing, 
body,  and  tail  aerodynamics.  Referring  to  Figure  2, 


-  sin  9 

sin  0  cos  0 
cos  0  cos  9 


F  A 


-D 

Y 

-L 


The  equations  determining  L  and  D  are  presented  later  in  this 
section.  Since  lateral  motion  is  not  addressed  in  this  paper,  sideforce 
(Y) ,  is  not  discussed. 

Tg  represents  the  engine  and  reaction  control  system  (RCS)  forces. 
For  simplicity  it  will  be  assumed  the  eg  and  all  the  RCS  and  main  engine 
nozzles  are  located  in  the  x-yg  plane.  All  RCS  thrust  vectors  will  be 
assumed  to  act  perpendicular  to  the  x-y  body  plane.  This  is  reasonable 
because  all  RCS  nozzles  are  8°  or  less  from  the  perpendicular.  The  main 
engine  nozzles  are  assumed  to  be  perpendicular  to  the  y  axis. 

D 
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Ib 


*“  — 

k 

r— '  — 

T 

T  cos  9. 

XB 

J 

T 

-T. 

3 

T 

-T  sin  9.  -  TQ  -  T.  -  T,  -  T, 

ZB 

j  6a  9f  *p  *s 

_  _ 

— 

The  RCS  thrusts  are  functions  of  the  stabilizer,  aileron  and  rudder 
positions  and  gross  thrust.  The  fore  and  aft  distribution  of  main 
nozzle  thrust  is  a  function  of  gross  thrust.  The  methods  for  computing 
the  thrusts  are  given  in  Chapter  2. 

Substituting  the  above  relations  into  equation  A-l  and  expanding. 


where  Ug, 


=  jj  [-  D  cos  a  cos  8  -  Y  cos  a  sin  8  +  L  sin  a  +  T  cos  9.]  -  sin  9 
w  J 

=  i  [-  D  sin  8  +  Y  cos  8  -  T^  ]  +  sin  0  cos  9  (A-2) 

»  ^  [■  D  sin  a  cos  8  -  Y  sin  a  sin  8  -  L  cos  a 

-  T  sin  9.  -  Tq  -  Ta  -  T.  -  T,  ]  +  cos  <j>  cos  9 

j  eA  9f  *p  *s 

and  w_  are  all  in  g's 

D 


Since  this  report  considers  only  flight  in  a  vertical  plane  and 
assumes  that  there  are  no  ambient  winds,  the  equation  for  v  drops  out 

D 

and  the  equations  for  u„  and  w  become 

D  D 


uD  =  w  ["  D  cos  a  +  L  sin  a  +  T  cos  0.]  -  sin  0 

D  n  J 


w_  =  7-7  [-  D  sin  a  -  L  cos  a  -  T  sin  0. 
B  W  j 


(A-3) 


Tq  -  Tq  ]  *  cos  9 
aA  °F 
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These  two  equations  are  the  basis  for  equations  F0(1)  and  F0(2)  in  the 
HARRIER. 

LIFT  AND  DRAG  EQUATIONS 

For  the  AV-8A  the  aerodynamic  lift  and  drag  functions  are  quite 

complicated  because  of  jet  interference  effects.  CT  is  not  only  a 

function  of  angle  of  attack  (a),  stabilizer  position  (6„),  flap  deflection 

n 

(<5p),  etc.  as  in  conventional  aircraft,  but  it  also  depends  on  the  nozzle 

deflection  (0^)  and  the  thrust  (T) .  Chapter  2  describes  how  and 

are  determined  based  on  available  data. 

1  2 

Given  C^,  C^.  and  C^,  and  q  A  y  pV^,  the  equations  for  aerodynamic 
lift,  side  force  and  drag  are  simply 
L  =  q  SCL 

Y  =  q  SCy  ( A - 4 ) 

D  =  q  SCD 

Again,  since  only  flight  in  the  vertical  plane  with  no  winds  is 
considered,  only  the  lift  and  drag  are  relevant  for  this  study. 

L  =  ^  SCL  (A- 5) 

D  =  q  SCD 

These  two  equations  are  the  basis  for  F0(3)  and  F0(4)  in  HARRIER. 

AXIS  TRANSFORMATIONS 

As  stated  in  Chapter  2,  inputs  are  defined  in  the  flight  path  axis 
system.  They  are  then  transformed  to  the  local  inertial  axis  system 
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for  comparison  to  the  aircraft  velocity  to  test  convergence. 


The  transformation  from  the  flight  path  axis  system  to  the  inertial 
axis  system  is  as  follows; 

X.  =  H1V  ( A  -  (> ) 

—I  v— v 

where 

*1  A  [Xj  Vj  tj] 

with  no  wind,  and  rotating  first  to  horizontal 


—  — 

— 

— 

—  — 

•s 

cos  Y 

cos 

c, 

-sin 

£ 

sin 

Y 

COS 

\' 

oo 

>1 

= 

cos  y 

sin 

c 

cos 

£ 

sin 

Y 

sin 

0 

:I 

-  sin 

■y 

0 

cos 

T 

0 

_ 

L_  — 

Where  is  the  aircraft  velocity  and  C  is  the  flight  path  azimuth 
angle  referenced  to  North.  Expanding  the  above; 


x,  =  V  cos  y  cos  £, 

I  00 

yT  =  V'  cos  y  sin  £ 

zT  =  -  V  sin  y 
I  00 

Assuming  £  to  be  zero,  can  be  omitted. 

x.  =  V  cos  y 

z,  =  -  V  sin  Y 
1  » 


(A- 7) 


(A- 81 


These  two  relations  express  the  commanded  inertial  velocity  components, 
VPC(l)  and  VPC(3),  in  the  equations  for  F0(5)  and  F0(6)  in  the  HARRIER. 


I2n 


. — 1 


mT"  mm . Mil 


The  remaining  components  of  equations  F0(5)  and  F0(6)  are  the 


elements  of  the  transformation  from  body  to  inertial  axes.  Given  , 

“*1 

the  components  of  the  wind  with  respect  to  the  earth. 


x  =  H*  V  +  V 
-I  nB  -B  -W, 


cos  0  cos  p  -  cos  <p  sin  t|/+  sin  d>  sin  iji+ 

sin  4>  sin  0  cos  <Ji  cos  <p  sin  0  cos  \p 

cos  0  sin  ^  cos  <j>  sin  ij;+  -  sin  $  cos  4>+ 

sin  <p  sin  0  sin  \p  cos  $>  sin  0  sin  \p 

-  sin  0  sin  <j>  cos  9  cos  $  cos  9 

Expanding  and  assuming  vD  *  <p  *  p  =  0  and  with  no  wind,  y  =  0  and 

o  1 


(A-9) 


“I 


!"•  — 

mm  _ 

UB 

V 

WXj 

VB 

+ 

V 

WVj 

WB 

— z\ 

M 

* 

>  _l 

— 

Xj  =  Ug  cos  9  +  Wg  sin  8 
ij  =  -  Ug  sin  0  +  Wg  cos  9 


(A-10) 


In  order  to  study  the  effects  of  acceleration  on  the  conversion 
corridor,  constant  accelerations  are  considered  inputs  in  the  flight 
path  axis  system.  These  were  transformed  to  the  body  axis  system  in 
which  the  aerodynamic  lift,  drag  and  moment  coefficients  are  available. 

Assuming  the  rate  of  change  of  wind  velocity  to  be  negligible. 


•  B  •  B  I  • 

V  =  H  V  =  H°  H  V 
-B  v  -v  I  v  -v 


(A— 11) 


Again,  <p  =  tp  =  £  =  vv  =  0,  so  for  our  problem  the  above  may  be  reduced  to 


Ug  =  uy  cos  (9-y)  -  wy  sin  (0-y) 


*B  =  uy  sin  (9-y)  +  wv  cos  (9-y) 


(A- 12 ) 
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However,  if  the  winds  are  also  assumed  to  be  zero 


a  =  0  -  y 

u„  =  u  cos  a  -  w  sin  a 
B  v  v 

Wg  =  Uy  sin  a  +  Wy  cos  a  (A-13) 

In  HARRIER,  uy  and  w^  are  referred  to  as  ACPC(l)  and  ACPC(3),  respectively, 
which  represent  the  inputs  APAR  and  APER. 


MOMENT  EQUATIONS 

In  general  for  a  rigid  body,  +  M^  =  ^  (IjWj)  in  an  inertial 

a  e 

coordinate  system,  where 

Mj  h  total  aerodynamic  moment 

a 

Mj  A  total  moment  due  to  engine  and  RCS  forces 

e 

IT  A  aircraft  inertia  tensor 

ojj  A  aircraft  rotational  velocity  vector 

In  quasi-equilbrium  flight  uij  =  0,  so  the  equations  reduce  to 

M  +  Mx  =  0 
a  e 

The  sum  of  the  moments  must  be  zero  irrespective  of  the  coordinate  system. 
Both  the  aerodynamic  and  the  engine  moments  are  given  in  body  axes. 


M 


B 

a 


+  M 

B 

e 


0 


M 

B 

a 


’l 

’l 

a 

e 

M 

M  A 

M 

a 

B  ” 

e 

N 

e 

N 

a_ 

e_ 

(A- 14) 


aerodynamic  moments 


engine  moments 


Assuming  the  moment  reference  center  is  located  along  the  Xg  axis,  the 
aerodynamic  moments  may  be  defined  as. 


La  =  q  Sb  C£ 

M  =  q  ScC  +  x  q  SC.  cos  a  +  x  q  SCn  sin  a 
a  n  m  mr  L  mr  ^  D 

N  =  q  SbC  +  x  q  SCV  cos  a 
a  n  n  mr  ^  Y 

The  components  used  in  computing  Cm  and  are  discussed  in  Chapter  2. 
C^,  Cy  and  are  not  discussed  in  this  paper. 

Again  referring  to  Figure  2, 


L  =  (T.  -  T ,  )  x 

e  <i>s  <t>p  <J> 


M  -  Ta  Xa  -Ta  xfl  +  T .  x._  -  T..  x.. 

6  0F  °F  0A  8A  jF  jF  jA  jA 


+  (T.  +  T.  )  x.fl 

4>p  <f>s  <J>9 


N  =  T,  x, 
e  ip  [f> 


The  modelling  of  the  RCS  and  engine  thrusts  is  described  in  Chapter  2, 
Substituting  the  above  into  equation  (A-14)  and  expanding. 


Rolling  0  =  q  SbC£  +  (T  -  T  )  x 

Moment 


Pitching  0  =  q  ScCm  ♦  ^qfX0F  -  TQa  xQA  -  T  x  (A-15) 

Moment 

-  T..  x..  +  (T,  +Ta  )  x.Q  +  x  q  SC.  cos  a  +  x  q  SCn  sin  a 
jA  jA  v  <pp  <PS  <pd  mr  ^  L  mr  D 

Yawing  0  =  q  SbCN  +  T  x  +  xmr  q  SCy  cos  6 

Moment  ‘  v  y 

As  previously  mentioned,  only  the  longitudinal  plane  is  considered 
in  this  paper  and  the  equation  for  pitching  moment  is  the  only  one 
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implemented  in  HARRIER.  In  addition  it  is  assumed  that  the  roll  and  yaw 
RCS  jets  are  not  being  used.  Therefore, 


Pitching  0  =  q  ScC  *  T  x  -  T  x  .  TjF  xjF 
Moment  F  F  A  A 

-  T..  x..  ♦  x  q  SC.  cos  a  +  x  a  SC..  sin  a 
jA  jA  mr  n  L  mr  M  D 


(A- 16) 


This  forms  the  basis  for  F0(7)  in  HARRIER. 


LI FT- LOSS  CURVES 

As  mentioned  in  Chapter  2,  some  discrepancy  was  found  in  the  data  in 
Reference  5  when  the  data  from  the  aforementioned  carpet  plots  for  C^  were 
compared  to  another  figure  in  the  same  Reference  (reproduced  here  as  Fig¬ 
ure  33)  which  shows  the  lift-loss  as  a  function  of  V./V  .  The  lift-loss 

j  OO 

presented  by  the  carpet  plot  data  in  Figures  3  through  6  can  be  compared 
to  that  in  Figure  33  in  the  following  way.  Assuming  the  lift- loss  is  in¬ 
variant  with  a  and  using  a  =  8°  as  reference, 

ACL  =  CL  -  CL 

wing-body  wing-body  °wing_body 

where 

C.  *  C,  at  given  0.  and  V./V^ 

Lwing-body  Lwing-body  J  ^ 

(from  Figures  3  through  6) 


Nondimens ionali zing  by  gross  thrust,  T; 


ACl  qS 

AL  _  wing-body 
T  T 

From  Reference  6,  T  =  (1/2) pV^.  =  (V.. /VJ 2qS_. 

AC,  qS  AC  S 

Therefore  ^  -  _ »in&-body  .  _ wing-boji^ 

(Vj/Vj'qSj  (Vj/VjV 

Using  the  above  formula,  the  lift- loss  presented  by  Figures  3  to  6  was  com¬ 
puted  and  superimposed  on  Figure  33.  These  data  obviously  do  not  agree.  Lift- 
loss  obtained  from  the  carpet  plots  is  substantially  greater  than  the  lift-loss 
graph  shows.  In  addition,  the  lift-loss  graph  shows  a  lift-loss  at  0^  =  0  which 
disagrees  with  the  observations  on  pages  12  and  13.  Furthermore,  as 

V./V  goes  to  zero,  the  lift-loss  curves  first  increase  and  then  decrease  in 
1  00 

a  range  where  no  substantiating  data  could  be  found.  The  reason  for  this 
might  be  that  the  lift- loss  should  go  to  zero  as  (hence  V^/VJ  goes  to 
zero.  This  would  be  reasonable  except  for  the  fact  that  at  this  point  both 
AL  and  T  tend  towards  zero  so  the  ordinate  of  this  figure,  AL/T,  is  actually 
an  ill-defined  point.  The  reason  for  these  discrepancies  could  not  be  de¬ 
termined. 

The  wing-body  aerodynamics  for  flaps -up  flight  were  taken  from 
Reference  5  and  are  presented  in  Figures  34  through  36.  These  data  were  ob¬ 
tained  with  =»  0  and  it  is  assumed  that  there  is  no  wing/ jet  interaction 
at  this  point. 
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APPENDIX  B 


METHOD  OF  SOLUTION 


In  earlier  work  in  this  field  Marinucci  (Ref.  10)  used  the  Bellman 
quasi -linearization  method  (Ref.  11)  for  a  similar  problem.  This  method 
is  an  extension  of  the  well  known  Newton-Raphson  algorithm.  Given  a  non¬ 
linear  function, 
y  =  f  (x) 

the  basic  Newton-Raphson  method  is  an  iterative  process  where  given  an 
initial  guess  for  y,  yQ  the  value  of  y  at  some  point  x^ew  may  be 
estimated  by, 


new  ~  ^old  dx 


df(x) 


x=x 


.  (x  -  x  .  ,) 
new  old 


old 


Bellman's  quasi-linearization  applies  this  to  a  system  of  equations,  for 


L  •  l  W 


^new  ~  ^old  3 


3f(x) 


x=x  ,  , 
—  -old 


(x  -  x  .  ,) 
-new  -old 


The  Jacobian  for  the  set  of  equations  used  in  this  problem  and 
typical  values  are  at  the  end  of  this  Appendix.  If,  as  in  our  case. 


y  ^  0  and  the  value  for  x  which  satisfies  this  is  desired,  then 

3f(x) 

...  -  x  .  .1  where  J  A 

-new 


0  =  y  ,  ,  +  J  (x  -  x  ,,) 
—  Hsld  —  -new  -old 


W  U  ft - 

—  -dx 


—new  *  ^old  ”  i  Xpid 

This  will  be  an  approximation  for  the  true  answer.  An  obvious 

condition  for  using  Bellmans'  method  is  that  J  be  non-singular.  Refs. 
10  and  11  describe  conditions  for  convergence  of  this  method,  the  most 


important  condition  being  the  convexity  of  the  functions  £(x) .  While 
convexity  of  the  functions  is  a  sufficient  criterion  for  convergence  it 
is  not  strictly  necessary.  For  less  dominant  terms  in  any  particular 


equation,  if  f.  (x.)  is  not  convex,  convergence  may  still  normally  be 
k  j 

obtained  by  setting  that  respective  term  in  the  Jacobian  to  zero. 

In  summary,  the  iterative  Bellman  quasi-linearization  is  used  as 


follows: 

1.) 

guess  a  reasonable  solution  x^^ 

2.) 

compute  the  Jacobian  J  1 

l£*2oid 

3.) 

solve  for  x  =  x  ,  ,  -  J_1y  ,  , 

-new  -old  —  ^old 

4.) 

test  convergence 

a. )  if  convergence  passes  -  STOP 

b. )  if  convergence  fails, 

make  x  .  ,  =  x  and  go  to  step  2. 

-old  -new  °  r 

Since  the  aerodynamic  and  control  system  functions  of  the  HARRIER 
are  highly  non-linear  most  of  the  functions  are  input  to  the  computer 
program  as  table  look-ups.  The  computer  SUBROUTINE  which  computes  the 
values  of  the  functions  also  computes  the  partial  derivatives  for  sub¬ 
stitution  into  the  Jacobian. 

With  the  equations  as  defined  in  Chapter  2,  the  Jacobian  is  shown 
on  the  last  pages  of  this  Appendix.  The  derivation  of  all  elements  in  the 
Jacobian  is  straightforward  except  for  the  partial  derivatives  with  respect 
to  T/W  since  the  aerodynamic  data  are  given  as  functions  of  V./V^,  and 


the  table  look-up  routine  therefore  gives  the  partial  derivative  with 


respect  to  V^/V^  ;  the  relation  between  V^/V^  and  T/W  must  be  known. 

T  =  (Vj/VJ2  q  S.  (as  defined  in  Ref.  6  ). 

The  derivative  with  respect  to  thrust  can  then  be  expressed  in  terms  of  the 
derivative  with  respect  to  V^/V^.  Given  the  partial  of  some  function  f^ 
with  respect  to  thrust. 


3f. 


3f. 


•  V.  2 

W  J  3(^1  qS. 

QO  J 


9f . 

l 


El  3(V 

w  W  J 


3f . 


However,  in  most  cases  the  table  look-up  computes  — — —  ,  therefore 


3  0 


3f.  3f . 

l  l 


V. 

30 

00 


3ff  5Si  v,  Vi  2 

V-  30 


3f . 

l 


qS.  V. 

—l  a  (_l) 

w  Vv  ; 


3fi  _  _W_ 
T  =  - 

4  ^ 


i  (li)-1 

'  2  V 


3f . 


V. 

2(^3 

00 


V. 

3(^3 


The  next  three  pages  show  the  elements  A(i,j)  of  the  Jacobian  in 
symbolic  form, 
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A 


* 


i 


A(l,l)  =  sin  a 

A(l,2)  =  -  cos  a 

A(l,3)  =  cos  0^ 


